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SIMPLIFIED METHOD FOR EXTRACTING MODEL PARAMETER SETS 
AND METHOD FOR STATISTICALLY SIMULATING INTEGRATED 
CIRCUIT USING THE SAME 

5 BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to a process of designing an integrated circuit, and 
more particularly, to a method for extracting a plurality of model parameters required for 
simulating the operation and the performance of a circuit design and a method for 
10 simulating the operation and the performance of the circuit. 

2. Description of the Related Art 

It is essential to determine whether integrated circuit designs will operate as 
desired, how the performance of the integrated circuit will be, and how the statistical 
distribution and deviation of the performances will be, and to feed back the simulation 

1 5 results into the design process. Such a simulation is performed using a widely known 
program, which is referred to as SPICE. The SPICE program requires information on 
the layout of devices such as the sizes of the circuit devices, the length of signal paths 
therebetween and various model parameters showing the characteristic performance of 
the devices and the overall interconnected circuit. Model parameters may have physical 

20 meaning such as the length xi and the width x w of a channel, the thickness t ox of a gate 
insulating layer, and the change AV t h in a threshold voltage according to the change in 
the doping concentration of a channel region. Other model parameters may indicate 
simple coefficients of a SPICE model equation. The process of extracting model 
parameters for a circuit is a complicated and time-consuming task, with many 

25 measurements and much trial and error. The general process of extracting model 
parameters will be described as follows, taking a transistor device as an example. 

Drain currents with respect to gate-source voltages Vgs and drain-source voltages 
V DS are measured and an I-V characteristic curve is prepared. The I-V characteristic 
curve can be measured, while changing the operating voltage Vcc, the substrate bias 

30 voltage Vbb, the length and the width of a channel, and device temperature. As a result, 
various I-V characteristic curves are created. 

In the I-V characteristic curves shown in FIG. 1, results of measuring drain 
current Id according to the change in the gate-source voltage Vgs> the drain-source 
voltage V D s, the substrate bias voltage Vbb and the operating voltage V C c> are shown. In 
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FIG. 1, 1-V characteristic curves with respect to predetermined length and width of a 
channel and a predetermined temperature are shown. Various I-V characteristic curves 
which are similar to the I-V characteristic curves shown in FIG. 1 are created as the 
length and the width of the channel and the temperature are changed. Because of all the 
5 variables, the number of times that a drain current value is measured can be several 

hundreds through several thousands. Typically, the transistor devices for measuring the 
current value according to the change in a voltage are referred to as a test element group 
(TEG) formed in an unused area, e.g., a scribe line between one chip and another chip, of 
a wafer on which mass-produced integrated circuit chips are formed. 

1 0 The process of extracting model parameters includes 1 ) substituting initial values 

into the model parameters for the I-V calculating equation of the SPICE program; 2) 
performing the needed calculations based on such initial values; and 3) changing the 
initial values until the calculation results converge on actual (measured) I-V 
characteristic curves. These steps are repeated until convergence is within a 

1 5 predetermined tolerance. Values obtained when the calculation results converge on the 
actually measured I-V characteristic curves become the model parameters. The model 
parameters of a previously known device having the characteristics closest to those of the 
currently designed device are typically used as initial values. The model parameters are 
preferably a set which converges on the measured I-V characteristic curves, while 

20 changes to the length and the width of the channel, the temperature during measurement, 
an operating voltage, and a substrate bias voltage concurrently are made. 

The model parameter set obtained as described above is substituted for the input 
variable to the SPICE model equation, and device and circuit performance 
characteristics, e.g. speed, are calculated. 

25 In order to determine the degree and deviation to which the performance of a 

circuit may vary according to the position of the circuit on a wafer or the position of a 
wafer in a wafer cassette, simulation is performed not on only one device but on plural 
devices variously positioned on one or more wafers (refer to US Patent No. US 
5,790,436). In this case, the above mentioned process of extracting the model parameter 

30 set must be repeated with respect to each device and wafer position to create a plurality 
of model parameter sets. Those of skill will appreciate that this is a complicated and 
time-consuming task. 

A method of extracting the model parameter set using electrical test (ET) data 
instead of I-V characteristic curves has been proposed. See James C. Chen et al "E-T 
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Based Statistical Modeling and Compact Statistical Circuit Simulation Methodologies", 
IEDM, 1996, pp. 635-638. That article describes a method of directly measuring 
parameters such as the change amount AL in the length of the channel, the change 
amount AW in the width of the channel, the thickness T 0 x of a gate insulating layer, and 

5 the drain-source resistance Ra s among the model parameters, and extracting the 

remaining parameters using ET data such as drain saturation current Idsat- However, it is 
more difficult to directly measure physical parameters than it is to measure the I-V 
characteristics. For example, it is very difficult to measure such parameters by a 
physical technique involving the use of a scanning electronic microscope. Accordingly, 

10 such parameters are usually measured by an indirect method by which the parameters are 
derived from electrical measurement. It is not possible to guarantee the reliability of 
such derived parameter values. 

Of interest to integrated circuit designers is what influence do changes in specific 
characteristics of an elementary device of an integrated circuit design have on the 

15 performance of the entire integrated circuit. Also of interest are which characteristics of 
the device most greatly affect the performance of the integrated circuit. However, 
conventional techniques for extracting model parameters and evaluating integrated 
circuit designs do not adequately meet such needs. 

20 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a method of simply extracting 
model parameter sets when two or more model parameter sets are required in order to 
statistically simulate the performance of a circuit design. 

It is another object of the present invention to provide a method of estimating 
25 sensitivity, i.e. determine the influence that changes in elementary device characteristics 
have on the performance of the entire circuit. 

Accordingly, to achieve the first object, there is provided a method of extracting 
n (n being an integer greater than or equal to 2) model parameter sets in order to 
statistically simulate the operation and performance of an integrated circuit design, 
30 wherein a first model parameter set is extracted by a conventional method and the 

remaining n-1 model parameter sets are extracted using not an I-V characteristic curve 
but instead the main characteristic data of an elementary device as a target function. 

In a method of extracting model parameter sets according to the present 
invention, n points are designated in one or more wafers and an elementary device of the 
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integrated circuit is formed at each of the n points. An I-V characteristic curve is created 
by measuring the I-V characteristic of the elementary device formed at one point among 
the n points and a first model parameter set is extracted which converges on the I-V 
characteristic curve within a predetermined error range. With respect to the second 
through nth model parameter sets, the second through nth main characteristic data values 
of the elementary device formed at the remaining second through nth points are 
measured and the second through nth model parameter sets are extracted which converge 
on the main characteristic data values within a predetermined error range. 

The elementary device may, for example, be a transistor and the main 
characteristic data may be one or more selected from the group consisting of a threshold 
voltage Vth, a drain saturation current I Ds at, a maximum conductance Gmm* x , a drain 
current W (when a gate-source voltage is OV), a junction capacitance Q, and an overlap 
capacitance C ov of the transistor. 

Also, the model parameters of the first model parameter set can be input as initial 
values for extracting the second through nth model parameter sets. 

According to the method of simulating the integrated circuit according to the 
present invention, it is possible to analyze the influence that changes in the main 
characteristic data value have on the performance of the integrated circuit. The model 
parameter sets are extracted using the main characteristic data as the target function. 

In a method of simulating the operation and the performance of a designed 
integrated circuit according to the present invention, an elementary device of the 
integrated circuit is formed at a predetermined point on a wafer. After creating an I-V 
characteristic curve by measuring the I-V characteristic of the elementary device, a first 
model parameter is extracted which converges on the I-V characteristic curve within a 
predetermined error range. A main characteristic data value of the elementary device is 
measured and the measured main characteristic data value is arbitrarily changed. A 
second model parameter set is then extracted which converges on the changed main 
characteristic data value within a predetermined error range. The first and second model 
parameter sets are input to a simulation program such as SPICE, the performances of the 
integrated circuit are simulated, and the simulated performances are compared with one 
another. Accordingly, the influence that the change in the main characteristic data value 
of the elementary device has on the performance of the integrated circuit design may be 
analyzed. 



PATENT APPLICATION 



4 



DO. NO. 9898-168 



Here, the first model parameter set can be extracted using not the I-V 
characteristic curve but the main characteristic data as the target function. 

According to the method of extracting n model parameter sets of the present 
invention, the simple-to-measure main characteristic data of the device (rather than the 
5 difficult-to-measure I-V characteristic curve) is used as the target function when those 
model parameter sets that follow the first model parameter set are extracted. Therefore, 
it is possible to significantly reduce the time spent on extracting the model parameter 
sets. 

Also, the invention makes it possible to easily analyze the influence that changes 
10 in the characteristics of the device have on the performance of the circuit. This is 

because the main characteristic data of the device are used to extract the model parameter 
sets. 



BRIEF DESCRIPTION OF THE DRAWINGS 
15 The above objects and advantages of the present invention will become more 

apparent by describing in detail a preferred embodiment thereof with reference to the 
attached drawings in which: 

FIG. 1 shows examples of I-V characteristic curves which are created by 
measuring current values according to the change in voltages in order to extract model 
20 parameter sets; 

FIG. 2 is a flowchart showing a method of extracting model parameter sets 
according to the present invention; 

FIG. 3 tabulates a 1st model parameter set and a kth model parameter set, 
extracted according to the present invention; 
25 FIG. 4 is a circuit diagram showing a ring oscillator as an exemplary integrated 

circuit design to be simulated; and 

FIG. 5 is a correlation diagram showing results of simulation of the ring oscillator 
using model parameter sets extracted by the conventional method and by a method 
according to the present invention. 

30 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
FIG. 2 is a flowchart showing the process of extracting n model parameter sets 
according to the present invention. FIG. 3 shows a 1st model parameter set and a kth 
model parameter set, extracted according to the method illustrated in FIG. 2. 
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Referring first to FIG. 2, a method of extracting n model parameter sets according 
to the present invention includes a extracting a first model parameter set, steps 100 
through 125, and a process of extracting a second through nth model parameter sets, 
steps 200 through 225. In the process of extracting the first model parameter set (steps 
5 100 through 125), a measured I-V characteristic curve is used as a target function as is 
conventional, but the process of extracting the second through nth model parameter sets 
(steps 200 through 225) constitutes a novel technique according to the present invention. 

N different points on one or more wafers are designated and an elementary device 
of a integrated circuit design such as a MOS transistor is formed at each designated 
10 point, in order to statistically simulate the operation and the performance of the 
integrated circuit. 

An I-V characteristic curve is provided by measuring the I-V characteristic of the 
device formed at an arbitrary point (a first point) among the n points (step 100). It is 
desirable to create various characteristic curves by measuring the I-V characteristic, 
15 while changing the length and the width of the channel of the device, the temperature 
during measurement, the operating voltage, and for the substrate voltage, thereby to 
improve the accuracy of the simulation. 

Next, initial values representing a first model parameter set are input (step 105). 
In this step, a previously known model parameter set of a device similar to the target 
20 elementary device design is selected, input or default values are provided by a SPICE 
program, and both are input as the initial values. 

An I-V characteristic value is calculated by substituting the initial values for the 
I-V calculating equation of the SPICE program (step 1 1 0). 

It is determined whether the I-V characteristic value calculated in step 110 
25 converged on the I-V characteristic curve measured in step 100, i.e. whether the initial 
values appropriately reproduce the measured I-V characteristic curve (step 115). It will 
be understood that the I-V characteristic value calculated in this step typically will not 
exactly coincide with the measured I-V characteristic curve. Accordingly, it is 
considered that the I-V characteristic value has converged on the I-V characteristic curve 
30 when the difference between the I-V characteristic value and the I-V characteristic curve 
is within a predetermined error range. 

Typically, the I-V characteristic value calculated by the original initial values 
input in the step 105 does not converge to the measured I-V characteristic curve. 
Accordingly, the initial values input in step 105 typically are changed (step 120). Then, 
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the I- V characteristic value is calculated again by substituting the changed initial values 
into the I-V calculating equation (step 110). In the process of changing the initial values 
(step 120), differences between the calculated I-V characteristic values and the measured 
I-V characteristic curves are evaluated. For example, it maybe that the calculated I-V 
5 characteristic value coincides with the measured I-V characteristic curve in the linear 
region of the I-V characteristic curve and the calculated I-V characteristic value is quite 
different from the measured I-V characteristic curve in the saturation region of the I-V 
characteristic curve. Or, conversely, it may be that the calculated I-V characteristic value 
coincides with the measured I-V characteristic curve in the saturation region of the I-V 

10 characteristic curve and the calculated I-V characteristic value is quite different from the 
measured I-V characteristic curve in the linear region of the I-V characteristic curve. 
Appropriate parameters (initial values) are changed according to the nature and extent of 
the difference. Preferably, the parameters are automatically changed and re-calculated 
by a known sub-program of the SPICE program (for which a plurality of optimization 

1 5 algorithms have been provided) and the Levenberg-Marquardt (L-M) algorithm using a 
least square method. (Refer to N. Arora, MOSFET Models for VLSI Circuit Simulation- 
Theory and Pranctice, Springer- Verlag, 1993, pp. 510-535). 

When the I-V characteristic value calculated in the step 110 converges on the 
measured I-V characteristic curve, the initial values which are finally input to the I-V 

20 calculating equation are extracted as a first model parameter set (step 125). In other 
words, a first model parameter set 300 shown in FIG. 3 is created. Each of the model 
parameter sets 300 and 400 shown in FIG. 3 is divided into two parts, that is, a part 
relating physical attributes 305 and 405 (the physical attributes of the parameters are 
described above) and a simple coefficient part 310 and 410 of the SPICE model equation 

25 (Refer to HSPICE User's Manual, Meta-Software Inc., 1996, pp. 16-103 through 16-106; 
and with respect to the detailed meanings of the illustrated coefficients). 

In the process of extracting model parameter sets with respect to the devices 
formed at the remaining n-1 points, the model parameter sets preferably are extracted by 
using several measured main characteristic data as a target function (rather than using the 

30 measured I-V characteristic curve). 

Namely, the main characteristic data values of the devices formed at a second 
point which is not the first point among n points are measured (step 200). ET data best 
representing the characteristic of the device, e.g. threshold voltage Vth, saturation current 
lDsat> a maximum conductance G mmzx , drain current when Vgs = 0V (Ioff), junction 
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capacitance Cj, and overlap capacitance C ov belong to the main characteristic data. These 
data can be obtained directly by a simple electrical measurement without resort to 
scanning electronic microscope or a complicated evaluation process. This distinguishes 
main characteristic data from, for example, the amount of change AL in the length of the 
5 channel, the amount of change AW in the width of the channel and the thickness T ox of a 
gate insulating layer. Further, the main characteristic data can be obtained by relatively 
few measurements. Therefore, it takes much less time and it is much easier to measure 
the main characteristic data than the I-V characteristic curve. 

The result of extracting the model parameters using only the main characteristic 

10 data which can be simply measured as the target function is not significantly different 
from the result of extracting the model parameters using the I-V characteristic as the 
target function. First of all, the model parameter sets of the devices formed at various 
points are extracted and the operation and the performance of the integrated circuits 
formed at various points are simulated in order to observe the deviation of the 

15 performances of the integrated circuits according to their positions in a wafer or a wafer 
cassette. In this case, the difference in the performances according to the positions in the 
wafer or the wafer cassette can be determined from the main characteristic data without 
resort to I-V characteristic curve measurements. This is because the difference in the I- 
V characteristic curves according to the positions is almost directly represented by the 

20 difference in the main characteristic data values of the device according to the positions 
(the first model parameter set having been extracted by using the measured I-V 
characteristic curve). 

The initial values which will become the second model parameter sets are input 
(step 205). In this step, the previously extracted first model parameter set (300 of FIG. 

25 3) is preferably directly input as the initial values. This is because the model parameter 
sets corresponding with the second through nth points are expected to be similar to the 
first model parameter set. 

The main characteristic data value is calculated by substituting the initial values 
for the calculating equation of the main characteristic data value of the SPICE program 

30 (step 210). 

It is determined whether the main characteristic data value calculated in step 210 
has converged on the main characteristic data value measured in step 200 (step 215). 
When the difference between the main characteristic data value calculated in step 210 
and the main characteristic data value measured in step 200 is within a predetermined 
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error range, it is concluded that the main characteristic data value calculated in step 210 
has converged on the main characteristic data value measured in step 200. 

When it is determined that the main characteristic data value calculated in step 
210 has not converged on the main characteristic data value measured in step 200, the 
5 initial values are changed (step 220). Then, the main characteristic data value is 

calculated again by substituting the changed initial values for the calculating equation of 
the main characteristic data value (step 210). In the step of changing the initial values 
(step 220), the above-mentioned optimization program can be used. Also, it is possible 
to change only parameters 305 having physical attributes without changing all the initial 

10 values (305 and 310 ofFIG. 3). This is because the parameters 305 having physical 
attributes among the model parameters change according to the conditions of actual 
processes and the remaining simple coefficient part 310 need not be changed according 
to the change in the positions in a wafer or wafer cassette. This is because the 
coefficients 310 have little relation to the change in the conditions of actual processes. 

15 When the main characteristic data value calculated in step 210 has converged on 

the measured main characteristic data value, the initial values finally input to the 
calculating equation of the main characteristic data value are extracted as the second 
model parameter set (400 of FIG. 3). As shown in FIG. 3, among model parameters of 
the second (k=2) model parameter set 400, model parameters which are different from 

20 those of the first model parameter set 300 are only the parameters 405 having physical 
attributes. This means that it is easier to create the target function and it is advantageous 
to the convergence of a solution to consider only those parameters (less than 10) having 
physical attributes according to the change in the processing conditions in the optimizing 
step (step 220). 

25 The third through nth model parameter sets are extracted by repeating steps 200 

through 225. 

Designating the first point for extracting the first model parameter set is arbitrary. 
However, after measuring the main characteristic data value of all the devices formed at 
n points, the device at the point where the main characteristic data value is closest to the 
30 specification of the integrated circuit design can be a reference sample (the first point). 

According to the method of extracting n model parameter sets according to the 
present invention, it is possible to significantly reduce the time spent extracting the 
model parameter sets. This is accomplished by extracting the model parameter sets by 
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measuring the I-V characteristic once and measuring only the main characteristic data 
(which can be easily measured the remaining (n-1) times), instead of by measuring the 
I-V characteristic (which takes a long time to measure and which is difficult to measure). 
EXPERIMENT EXAMPLE 1 
5 The time delay through an 1 1 -stage ring oscillator formed of CMOS transistors as 

shown in FIG. 4 is simulated in order to check the accuracy of the simulation of the 
method according to the present invention. Transistor devices are formed at ten points 
across three wafers. Ten model parameter sets are extracted by the conventional method, 
in which the I-V characteristic curves are provided with respect to all the devices formed 

10 at the ten points. Other ten model parameter sets are extracted by the method according 
to the present invention, in which one model parameter set is extracted by measuring the 
I-V characteristic with respect to only the reference sample formed at one point and the 
remaining nine model parameter sets are extracted by measuring only a threshold voltage 
and a drain saturation current as the main characteristic data. Time delay is simulated by 

1 5 inputting the respective 10 model parameter sets extracted by the conventional method 
and the method according to the present invention. 

The results of simulating time delay with the model parameter sets extracted by 
the conventional method and by the method according to the present invention, with 
respect to four combinations of operating voltage Vcc and temperature, are shown in 

20 Table 1 . The mean error of the results of the time delay is 3. 1 5%, which is excellent. 



Table 1 



Simulation conditions 


Mean error (%) 


Vcc (V) 


Temperature (°C) 




2.5 


85 


3.3 


2.5 


110 


3.1 


3.0 


85 


3.5 


3.0 


110 


2.7 


Mean 


3.15 
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FIG. 5 is a graph showing the results of simulating the time delay with respect to 
the model parameter sets extracted by the conventional method and the method 
according to the present invention when V cc is 2.5V and the temperature is 85°C (the first 
entry in Table 1). Referring to FIG. 5, it may be seen that the result of simulating the 
5 time delay by the method according to the present invention closely coincides with the 
result of simulating the time delay by the conventional method. Namely, time delay 5 1 0 
simulated with respect to the reference sample using the first model parameter set 
extracted using the I-V characteristic curve as the target function and remaining nine 
simulated time delays are distributed closely along a diagonal straight line 520. Thus, 

1 0 the graph illustrates that the result of simulating the time delay by the method according 
to the present invention closely coincides with the result of simulating the time delay by 
the conventional method. The above assumption — that the difference in the main 
characteristic data values of the devices according to the positions on the wafer directly 
represents the difference in the I-V characteristic curves — is verified. 

1 5 Also, according to the method of extracting the model parameter sets of the 

present invention, it is possible to analyze the sensitivity, i.e. the influence that the ET 
data (the main characteristic data) has on the performance of the circuit. Namely, 
according to the present invention, it is possible to estimate the change in performance of 
the integrated design circuit according to the change in device characteristics without 

20 manufacturing the device whose characteristics (ET data) are changed. Further, it is 

possible to analyze the characteristic of the device that has the greatest influence on the 
performance of the integrated circuit design by comparing a) the result of simulating 
performance with respect to the model parameter set extracted using the I-V 
characteristic curve or the ET data as the target function with b) the result of simulating 

25 performance with respect to the model parameter set extracted using changed ET data as 
the target function. The following experiment examples show an example of such 
sensitivity analyses. 
EXPERIMENT EXAMPLE 2 

After extracting the model parameter set with respect to one reference sample and 

30 simulating the time delay of a 23-stage ring oscillator using the model parameter set, 

corresponding time delays using the model parameter sets extracted while changing the 
ET data are simulated, and the two are compared with each other. The simulated time 
delay of the reference sample was 3.28 nsec. The simulated time delays and the change 
ratios when the ET data is changed are shown in Table 2. In Table 2, C ja , Gj p , C ov , fosat, 
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and V t h denote a junction capacitance of source/drain regions to a channel, junction 
capacitance of the source/drain regions to a substrate, overlap capacitance, drain 
saturation current, and threshold voltage, respectively. For example, I Ds at-10% means 
that Idsat is reduced by 10% and the remaining characteristics are maintained. 
5 Table 2 





Changed ET data 


Time delay (nsec) 


Change ratio of time 
Uciay \/o) 


NMOS 


O I 1 CO/ 

Uj a + 1 J7o 




yj.D 




1 oo 


o a 


C ov + 15% 


3.32 


1.3 


IDsat" 10% 


3.50 


o.o 


v, +i n% 

V th ' iv/o 




0 1 

V/. 1 


PMOS 


Cja+15% 


3.33 


1.6 


Qp + 15% 


3.30 


0.7 


C ov +15% 


3.37 


2.9 


lDsat+15% 


3.62 


10.4 


V th -10% 


3.23 


0.2 



It is noted from Table 2 that the factor that has the greatest influence on the time 
delay of the ring oscillator is I Ds at, whether NMOS or PMOS transistors are used in the 
invention. 

1 0 EXPERIMENT EXAMPLE 3 

With respect to an asynchronous DRAM, after extracting the model parameter set 
of one reference sample and simulating tRAc (access time from when an address is 
designated to when data is read) using the model parameter set, tRAc is simulated using 
the model parameter sets extracted while changing the ET data, and the two are 

15 compared with each other. tRAc is 40 nsec when V C c is 3.0V and the temperature is 75°C 
on the average. The influence (the sensitivity) that the change in each ET data has on 
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tRAC is shown in Table 3. In Table 3, I Dp , fern, R pp , and C par a denote I Dsat of the PMOS 
transistor, I Dsa t of the NMOS transistor, the resistance of an upper electrode line formed 
of polycrystalline silicon, and the parasitic capacitance of metal wiring, respectively. 

Table 3 



ET data 


Influence that change in each ET data has 
on tRAc (sensitivity) (%) 


Idp 


36 


iDn 


21 


Rpp 


18 


C pa ra 


6 


Others 


19 



5 

It is noted from Table 3 that the factor that has the greatest influence on tRAc is 
the change in the saturation current of the PMOS transistor. Accordingly, a DRAM 
designer can estimate the influence that the changes in the characteristics of a DRAM 
have on the performance of an entire DRAM circuit without actually manufacturing a 

1 0 DRAM whose characteristics are changed. 

According to the present invention, in the case of extracting n model parameter 
sets in order to statistically simulate the performance of the integrated circuit, it is 
possible to significantly reduce the time spent on extracting the model parameters by 
extracting only one model parameter set using the I-V characteristic curve created by 

1 5 directly measuring the I-V characteristic of the device as a target function and extracting 
the remaining (n-1) model parameter sets, using the main characteristic data (ET data) 
that can be easily measured, e.g. the threshold voltage or saturation current of the device 
as the target function. 

According to the present invention, since the main characteristic data (ET data) of 

20 the device is used as the target function, it is possible to easily extract the model 

parameter set when the characteristics of the device are changed. Therefore, the designer 
can simply estimate the influence (sensitivity) that the changes in the characteristics of 
the device have on the performance of the integrated circuit. 
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Having illustrated and described the principles of my invention in a preferred 
embodiment thereof, it should be readily apparent to those skilled in the art that the 
invention can be modified in arrangement and detail without departing from such 
principles. I claim all modifications coming within the spirit and scope of the 
accompanying claims. 
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WHAT IS CLAIMED IS: 



1 . A method of extracting n wherein n is a natural number greater than or 
equal to 2, model parameter sets in order to statistically simulate the operation and the 
5 performance of an integrated circuit design, the method comprising the steps of: 

designating n points in one or more wafer and forming an elementary device of 
the integrated circuit at each of the n points; 

creating an I-V characteristic curve by measuring the I-V characteristic of the 
elementary device formed at one point among the n points; 
1 0 extracting a first model parameter set which converges on the I-V characteristic 

curve within a predetermined error range; 

measuring a main characteristic data value of the elementary device formed at a 
kth, wherein k is a natural number from 2 to n, point among the n points; and 

extracting a kth model parameter set which converges on the main characteristic 
1 5 data value within a predetermined error range, 

wherein n model parameter sets are extracted by repeating the steps of measuring 
the main characteristic data value until k is equal to n, increasing k by one, and 
extracting the kth model parameter set. 

20 2. The method of claim 1 , wherein the elementary device is a transistor and 

the main characteristic data are one or more selected from the group consisting of a 
threshold voltage V th , a drain saturation current I Ds at> a maximum conductance Gmmax, a 
drain current I off when a gate-source voltage is 0V, a junction capacitance C j? and an 
overlap capacitance C ov of the transistor. 

25 

3 . The method of claim 2, wherein the main characteristic data is the 
threshold voltage V th and the drain saturation current I Dsat of the transistor. 
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4. The method of claim 1 , wherein the step of extracting the kth model 
parameter set comprises the steps of: 

inputting the first model parameter set as initial values of model parameters 
which will become the kth model parameter set; 
5 calculating a main characteristic data value by substituting the initial values for 

the calculating equation of the main characteristic data value; 

determining whether the calculated main characteristic data value converges on 
the measured main characteristic data value within a predetermined error range; 

extracting the initial values as the kth model parameter set when it is determined 
10 that the calculated main characteristic data value converges on the measured main 
characteristic data value; and 

changing the initial values when the difference between the calculated main 
characteristic data value and the measured main characteristic data value deviates from 
the predetermined error range and repeating the steps of calculating the main 
1 5 characteristic data value and determining whether the calculated main characteristic data 
value converges on the measured main characteristic data value within a predetermined 
error range. 

5. The method of claim 4, wherein the initial values are changed only in 
20 model parameters having physical attributes. 

6. The method of claim 1, further comprising the step of arbitrarily changing 
the measured main characteristic data value of the elementary device, formed at an m, 
wherein m is an arbitrary natural number from 1 to n, point and extracting an mth model 

25 parameter set which converges on the arbitrarily changed main characteristic data value 
within a predetermined error range, in order to analyze the influence that the change in 
the main characteristic data value of the elementary device has on the performance of the 
integrated circuit design. 
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7. A method of simulating the operation and the performance of an 
integrated circuit design, the method comprising the steps of: 

forming an elementary device of the integrated circuit at a predetermined point 
on a wafer; 

5 creating an I-V characteristic curve by measuring the I- V characteristic of the 

elementary device; 

extracting a first model parameter which converges on the I-V characteristic 
curve within a predetermined error range; 

measuring a main characteristic data value of the elementary device; 
10 arbitrarily changing the main characteristic data value and extracting a second 

model parameter set which converges on the changed main characteristic data value 
within a predetermined error range; and 

inputting the first and second model parameter sets to a simulation program, 
simulating the performances of the integrated circuit based on the first and second model 
15 parameter sets, and comparing the simulated performances with each other, 

thereby analyzing the influence that the change in the main characteristic data 
value of the elementary device has on the performance of the integrated circuit design. 

8. The method of claim 7, wherein the elementary device is a transistor, the 
20 main characteristic data are one or more selected from the group consisting of a 

threshold voltage V t h, a drain saturation current Iosat, a maximum conductance Gm max> & 
drain current W when a gate-source voltage is 0V, a junction capacitance Q, and an 
overlap capacitance C ov of the transistor, each influence that the change in the main 
characteristic data value has on the performance of the integrated circuit design being 
25 analyzed by individually changing one value among the main characteristic data or 

combining two or more values with each other to produce change or combination results 
for purposes of extracting the second model parameter set. 

9. The method of claim 7, where the simulation program is SPICE. 
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10. A method of simulating the operation and the performance of an integrated 
circuit design, comprising the steps of: 

forming an elementary device of the integrated circuit at a predetermined point 

on a wafer; 

measuring a main characteristic data value of the elementary device; 

extracting a first model parameter set which converges on the main characteristic 
data value within a predetermined error range; 

arbitrarily changing the main characteristic data value and extracting a second 
model parameter set which converges on the changed main characteristic data value 
within a predetermined error range; and 

inputting the first and second model parameter sets to a simulation program, 
simulating the performances of the integrated circuit, and comparing the simulated 
performances with each other, 

thereby the influence that the change in the main characteristic data value of the 
elementary device has on the performance of the integrated circuit is analyzed. 

1 1 . The method of claim 10, wherein the elementary device is a transistor, the 
main characteristic data are one or more selected from the group consisting of a 
threshold voltage V th , a drain saturation current I Dsat , a maximum conductance G m max, a 
drain current I of f when a gate-source voltage is 0V, a junction capacitance Q, and an 
overlap capacitance C ov of the transistor, each influence that the change in the main 
characteristic data value has on the performance of the integrated circuit being analyzed 
by individually changing one value among the main characteristic data or combining two 
or more values with each other to produce change or combination results for purposes of 
extracting the second model parameter set. 

12. The method of claim 10, wherein the simulation program is SPICE. 
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SIMPLIFIED METHOD FOR EXTRACTING MODEL PARAMETER SETS 
AND METHOD FOR STATISTICALLY SIMULATING INTEGRATED 
CIRCUIT USING THE SAME 



5 ABSTRACT OF THE DISCLOSURE 

A method of simulating the performance of an integrated circuit design is 
provided. In the case of extracting n model parameter sets in order to statistically 
simulate the performance of the integrated circuit, it is possible to significantly reduce 
the time spent on extracting the model parameters by extracting only one model 

10 parameter set using the I-V characteristic curve created by directly measuring the I-V 

characteristic of the device as a target function and extracting the remaining (n-1) model 
parameter sets using the main characteristic data (ET data) which can be easily measured 
such as the threshold voltage or the saturation current of the device as the target function. 
Since the main characteristic data (ET data) of the device is used as the target function, it 

1 5 is possible to easily extract the model parameter set when the characteristics of the 
device are changed. Therefore, the designer can simply estimate the influence 
(sensitivity) that the changes in the characteristics of the device have on the performance 
of the integrated circuit design. 
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